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Recent  modeling  attempts  start  with  very  small  "point^rainfall  distribution 
functions  which  are  transformed  into  specific  attenuation  data.  The^point^ 
rainfall  rates  are  based  upon  a large  history  of  data  which  does  not  indicate 
the  spatial  variation  of  rain.  There  are,  however,  indications  that  outage  time 
on  1 ine-of-sight  communication  links  can  be  estimated  from  distributions  of 
L"'point**4'rainfall  rates, derived  from-U-.Sv  father  Service  information*. 

Due  to  the  random  nature  of  rainfall,  and  due  to  the  time  dependence  of  such 
information,  a logical  approach  to  forecasting  this  phenomenon  and  interpreting 
the  results  with  respect  to  systems  performance  seems  to  lie  within  the  realm  of 
non-stationary  time-series  modeling.  This  report  presents  an  attempt  to  develop 
statistical  models  which  can  be  used  to  forecast  in  near-real -time  and  to 
characterize  the  underlying  stochastic  processes  of  short-term  rainfall  density 
information.  It  should  be  mentioned  that  the  analysis  and  modeling  are  directed 
towards  the  accurate  characterization  of  precipitation  density  and  not  the 
probability  of  precipitation  occurrence. 

Due  to  the  wavelengths  involved  in  Line-of-Sight  Communication  links,  i.e., 
5-30  GHz,  the  size  of  raindrops  have  a definite  dispersive  and  absorptive  effect| 
on  propagated  electromagnetic  energy.  The  importance,  therefore,  of  this  work 
is  self-evident. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEPWlon  *m, rmd) 


CONTENTS 


! 


t 


t 

r- 
* t 

« 

r 

i 

» 

• « 

f 


p t 


r 

i 


i 


PAGE 


TIME  SERIES  MODELING  OF  CLIMATOLOGICAL  INFORMATION 

1.  INTRODUCTION 1 

2.  CONCEPTS  IN  TIME-SERIES 4 

2.1  Stationary  and  Non-Stationary  Time-Series 5 

2.2  Parametric  Time-Series  Models 7 

2.2.1  Selecting  the  Best  Model 8 

2.2.2  Estimation  of  Parameters, 8 

2.2.3  Checking  the  Fit 10 

2.2.4  Forecasting 11 

3.  IDENTIFICATION  OF  CLIMATOLOGICAL  DATA 11 

3.1  Fitting  the  Models 22 

3.2  Checking  the  Fit 29 

3.3  Proposed  Approach  vs.  Akaike's  Approach 29 

3.4  Summary  and  Conclusions 33 

BIBLIOGRAPHY 35 

FIGURES 

3.1  Measured  Daily  Precipitation  for  Greenwood  Lakes,  N.  J., 

January  through  December,  1974 12 

3.2  Measured  Daily  Precipitation  for  Long  Branch,  N.  J., 

January  through  December,  1974 13 

3.3  Sample  Autocorrelation  of  the  Daily  Precipitation  for 

Greenwood  Lakes,  N.  J.,  1974 15 

3.4  Sample  Autocorrelation  of  the  Daily  Precipitation  for 

Long  Branch,  N.  J.,  1974 16 

3.5  Sample  Autocorrelation  of  the  First  Difference  Data  of  the 

Daily  Precipitation  for  Greenwood  Lakes,  N.  J.,  1974 17 


i 


CONTENTS 


PAGE 


FIGURES  (cont'd) 

3.6  Sample  Autocorrelation  of  the  First  Difference  Data  of  the 

Daily  Precipitation  for  Long  Branch,  N.  J.,  1974 18 

3.7  Sample  Autocorrelation  of  the  Second  Difference  Data  of  the 

Daily  Precipitation  for  Greenwood  Lakes,  N.  J.,  1974 19 

3.8  Sample  Autocorrelation  of  the  Second  Difference  Data  of  the 

Daily  Precipitation  for  Long  Branch,  N.  J.,  1974 20 

3.9  Sample  Autocorrelation  of  the  Third  Difference  Data  of  the 

Daily  Precipitation  for  Long  Branch,  N.  J.,  1974 21 

3.10  Model  Order  vs.  Residual  Variance  for  the  Daily  Precipitation, 

1974,  for  Greenwood  Lakes,  N.  J 23 

3.11  Model  Order  vs.  Residual  Variance  for  the  Daily  Precipitation, 

1974,  for  Long  Branch,  N.  J 24 

3.12  Simulated  Precipitation  Density  Series  Using  the  Mixed  (ARMA) 


Model  vs.  the  Observed  Series  for  Greenwood  Lakes,  N.  J.,  1974....  26 


3.13  Simulated  Precipitation  Density  Series  Using  the  Mixed  (ARMA) 

Model  vs.  the  Observed  Series  for  Long  Branch,  N.  J.,  1974 27 

TABLES 

3.1  Kendall's  Tau  Statistics  for  Trend,  z0  = +1.645 14 

3.2  Forecasted  Values  of  Rainfall  Density  for  Greenwood  Lakes,  N.  J., 
(1974)  at  Origin  t = 142,  and  Updating  Under  the  Assumption  that 

x^43  Becomes  Available 28 

3.3  Forecasted  Values  of  Rainfall  Density  for  Long  Brand..  N.  J., 

(1974)  at  Origin  t = 144,  and  Updating  Under  the  Assumption  that 

X145  Becomes  Available 28 

3.4  Sample  Autocorrelation  of  the  Residuals,  rzz(k),  for  the 

Simulated  Precipitation  Density,  Greenwood  Lakes,  N.  J.,  1974 30 


3.5  Sample  Autocorrelation  of  the  Residuals,  rzz(k),  for  the 
Simulated  Precipitation  Density,  Long  Branch,  N.  J.,  1974. 

3.6  Best  Parameter  Order  for  the  Minimum  Residual  Variance  and 

Aka ike's  Criteria 


ii 


■ ■»  : — 


TIME-SERIES  MODELING  OF  CLIMATOLOGICAL  INFORMATION 


1.  INTRODUCTION 

It  is  known  that  communications  systems  operating  in  the  millimeter-wave 
range  may  be  critically  dependent  on  meteorological  factors.  Specifically, 
atmospheric  phenomena  such  as  hard,  steady  rainfall,  and  thunderstorms  may 
cause  serious  outages.  In  this  chapter,  the  phenomenon  of  rainfall  will  be 
addressed,  where  the  accuracy  of  weather  forecasts  will  not  be  considered. 
Simply,  the  interest  here  is  not  whether  or  not  there  will  be  precipitation, 
but  rather,  from  the  communicators'  point  of  view,  what  density  of  precipita- 
tion can  be  expected  if  it  does  rain.  Such  information  becomes  an  important 
tool  in  determining  the  need  for  back-up  modes  of  communications. 

Work  done  by  Lin,  [l  ],  and  Chen,  [2  ],  indicate  techniques  that  allow 
the  basis  of  long-term  information  for  empirical  conversion  of  rainfall  rate 
x ( t ) for  a given  integration  time  T on  the  basis  of  60-m.inute  rainfall  data. 
Their  technique  provides  a procedure  of  obtaining  a fairly  small  (5-minute) 
rain-rate  distribution  from  NOAA  climatological  data  on  short-time  rates  of 
rainfall  (i.e.,  thunderstorm  activity).  Obviously,  there  is  a great  stress 
on  the  "short-interval"  point-rainfall  distribution.  Rice  and  Holmberg, 

[3],  have  provided  a widely  accepted  model  utilizing  the  point-rainfall 
distribution  in  obtaining  the  "fraction  of  time  during  which  t minute  average 
rainfall  rates  exceed  any  given  value."  This,  in  turn,  is  used  to  project 
adverse  effects  of  rain  on  microwave  radio  links.  Stated  simply,  they 
postulate: 

rainfall  * mode  1 rain  + mode  2 rain  . 

That  is,  rainfall  density  is  composed  of  mode  1,  an  Individual  exponential 
mode  which  corresponds  to  a physical  analysis  of  thunderstorms,  and  that  of 
mode  2,  corresponding  to  all  other  rain.  The  number  of  hours  of  rainy  "t 
minute  periods"  for  which  an  average  surface-point  (station)  rainfall  rate 
is  exceeded  is: 

Tt(R)  = Tltqlt(R)  + T2tq2t^  hours  ’ (l*1) 
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where  Tlt  and  T2t  are  the  average  annual  rainfall  totals  for  modes  1 and  2, 
respectively,  and 

qu(R)  = exp  (-R/_  ), 

R1t 

the  time  that  a rate  R is  exceeded  by  mode  1 rain, 

q2t(R)  * 0.35  exp(-. 453074  R/  ) + 0.65  exp(-2. 857143  R/_  ) , 

*2t  R2t 

the  time  that  a rate  R is  exceeded  by  mode  2 rain, 

ftjt  = Mi^Tn  nOT^hr  » 
the  average  annual  mode  1 rainfall  rate, 

R?t  = M2^Tlt  tnni^hr  * 

the  average  annual  mode  2 rainfall  rate;  and  M2  are  mode  1 and  2 annual 
rainfall  totals,  respectively. 

Dutton  et  al , [4],  used  this  methodology  with  modifications  to  predict 
communications  link  performance  in  the  European  environment.  Among  their 
modifications,  a more  sophisticated  approach  in  estimating  the  parameters, 

T1t  T2t’  Rlt’  and  R2t’  usin9  multiple  linear  regression  was  proposed.  In 
both  of  these  aforementioned  works,  extensive  use  is  made  of  contour  maps 
of  annual  precipitation  with  appropriate  interpolation  to  estimate  parameters. 
Temporal  (year  to  year)  and  seasonal  variations  were  assessed  from  monthly 
precipitation  data.  Similar  work  by  Bartow,  [5],  is  also  based  on  long-term 
statistics  and  the  assumption  of  a constant  rain  rate.  These  aforementioned 
works  are  widely  used  to  predict  the  effects  of  weather  on  communications 
links  operating  in  the  8-30  MHz  band  and  basically  yield  approximations  of 
sizable  outage  times  (i.e.,  hours/year)  due  to  rain.  Sufficient  data  is  not 
yet  available  to  validate  the  prediction  error  [4  ] for  specified  confidence 
intervals. 

Since  the  available  data  is  random  and  time  dependent,  Jones,  [6  ], 

[7],  suggests  the  fitting  of  autoregressive  (AR)  models  to  the  given  time 
series.  The  critical  question,  of  course,  is  identification  of  the  order  of 
the  autoregressive  process.  In  addition,  one  is  always  concerned  about  the 
choice  of  the  autoregressive  model  over  the  other  two  commonly  used  models. 


namely,  the  moving  averages,  and  mixture  of  autoregressive  and  moving 
averages  models.  Jones  seems  to  ignore  these  models  in  his  research  of  the 
subject  area.  He  implies  in  a later  paper,  [8],  that  the  modified  Akaike’s 
FPE  criteria  is  useful  for  this  purpose.  However,  one  cannot  be  certain  as 
to  which  of  the  three  difference  equations  will  best  characterize  the  given 
rainfall  information.  Using  the  asymptotically  unbiased  estimate  of  mean- 
square  error: 


n 

where  S_  = Z ej,  n is  the  number  of  observations,  and  p is  the  number  of 
p t-1  1 

autoregressive  parameters,  he  obtains: 

FPEp  . ^ {1  ♦ (1  ♦P/n)>Ux*  . (1.3) 

This  criterion  is  computed  for  each  order,  p.  The  order  that  corresponds  to 
min(FPEp)  is  equivalent  to  the  minimum  variance  and  autocorrelation  criteria 
commonly  used.  Again,  no  mention  is  made  of  the  utility  of  the  integrated 
moving  averages  (IMA),  and  mixed  models  in  modeling  this  type  of  climatolog- 
ical information. 

In  view  of  these  shortcomings,  it  should  also  be  noted  that  the  time- 
series  approach  is  considered  here,  not  only  because  the  data  is  time 
dependent,  but  because  weather  observations  are  usually  made  some  distance 
away  from  the  area  of  use  and  may  not  be  fully  representative  of  conditions 
at  the  point  of  use,  and  there  is  usually  a time  lag  between  observations  and 
use.  These  factors,  along  with  random  occurrence,  contribute  to  the 
non-stationarity  of  the  underlying  process  in  that  the  values  of  the  climato- 
logical data  may  change  during  the  time  lag.  If  the  non-stationarities  are 
not  properly  approached,  then  meaningful  results  of  analysis  would  be  impos- 
sible to  obtain.  In  this  chapter,  the  time-series  methodology  will  be  shown 
for  daily  rainfall  information  for  Greenwood  Lakes,  New  Jersey,  and 
Long  Branch,  New  Jersey,  taken  from  the  1974  New  Jersey  Climatological  Data, 
National  Oceanic  and  Atmospheric  Administration,  U.  S.  Department  of 
Commerce,  [9  ].  These  data  are  the  station  accumulations  taken  from 
recording  rain  gauges.  Inference  will  be  made  as  to  the  applicability  of 
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associated  probability  distribution.  Thus  we  may  describe  the  behavior  of  a 
time-series  at  all  instances  by  an  ordered  set  of  random  variables  and  the 
associated  probability  distributions,  denoted  by  {Xt } and  fxt  , t = 0,  +1, 

+2,  •••.  Such  an  ordered  set  of  random  variables  is  called  a stochastic 
process.  Thus,  an  observed  time-series  xt  can  be  considered  as  one  realiza- 
tion of  an  infinite  ensemble  of  functions  which  may  have  been  generated  by  a 
stochastic  process.  A stochastic  process  is  said  to  be  strictly  stationary 
if  the  joint  distribution  of  any  set  of  observations  is  unaffected  by  shift- 
ing all  times  of  observations  ahead  or  backward  by  any  integer  amount  k.  A 
stationary  stochastic  process  may  be  described  in  terms  of  its  mean  u which 
is  estimated  by: 

x - n ^LjXt  > (2.1) 

its  variance  a2  which  is  estimated  by: 

Sx  = ; " , (xt  - x)2.  (2.2) 

x n t=l  L 

its  sample  autocovariance  function,  which  measures  the  extent  to  which  two 
random  variables  are  linearly  independent: 

cxx(k)  = |-^xt  " x)Ut+k  " x)  k = 0,  1,  •••,  n - 1,  (2.3) 

and  the  sample  autocorrelation  function,  which  acts  like  a correlation  co- 
efficient: 

rxx(k)  = cxx(k)/cxx(0)  k = 0,  1,  •••,  n - 1.  (2.4) 

2.1.  Stationary  and  Non-Stationary  Time  Series 
A stationary  time-series  is  one  which  is  in  statistical  equilibrium  in 
the  sense  that  its  properties  do  not  change  with  respect  to  time,  whereas  a 
non-stationary  time-series  is  such  that  its  properties  change  with  time. 

Time-series  occurring  in  practice  are  usually  non-stationary  in  nature  and 
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can  be  divided  into  three  classes: 


(i)  Those  which  exhibit  stationary  properties  over  a long  period  of 


time. 


(ii)  Those  which  are  approximately  stationary  over  very  short  periods 
of  time. 

(iii)  Those  which  exhibit  non-stationary  properties,  that  is,  their 
visual  properties  change  continuously  with  time. 

At  present  there  exist  techniques  to  analyze  stationary  time-series, 
but  the  techniques  available  for  the  analysis  of  non-stationary  time-series 
are  inadequate  and  do  not  lend  themselves  to  meaningful  interpretations  of 
physical  problems.  However,  one  can  adjust  non-stationary  time-series  so 
that  the  existing  techniques  of  stationary  time-series  analysis  can  be 
applied.  The  adjustment  is  accomplished  by  applying  a proper  filter  to  the 
observed  non-stationary  time-series  to  remove  the  non-stationary  components. 

The  selection  of  a proper  filter  is  accomplished  through  a search  for  a 
mathematical  function  which  will  transform  a non-stationary  series  into  a 
stationary  series.  One  of  the  most  used  and  most  efficient  methods  of 
filtering  is  through  the  application  of  a difference  equation,  [13],  [14]. 


A first-order  difference  equation  is  defined  by: 
n = xt  - xt-l» 


(2.5) 


where  xt  is  the  observed  non-stationary  series  and  yt  is  the  first-difference 
series.  Similarly,  a second-order  difference  equation  is  defined  by: 

wt  = xt  " 2xt-l  + xt-2’  <2-6) 

and  so  on.  A first  or  second-order  difference  equation  will  usually  be 


sufficient  to  transform  most  practically  occurring  non-stationary  time-series, 

P31- 


To  identify  whether  the  observed  series  exhibits  stationary  or  non- 
stationary properties  one  can  use  certain  data-analysis  tools.  In  addition 
to  graphical  representation  of  the  observed  series  the  sample  autocorrelation 
function  of  the  observed  series  and  a trend  test  applied  to  the  observed 
series  are  important.  For  the  observed  series  and  its  first  and  second 
differences,  the  sample  autocorrelation  functions  (2.4)  are  computed  and  a 
trend  test,  Kendall's  tau,  £llj  , is  performed.  (The  sample  autocorrelation 
function  of  a stationary  series  has  the  property  that  is  dampens  out  fairly 
rapidly,  that  is,  it  approaches  zero;  also,  a stationary  series  will  be  such 
that  it  contains  no  trend).  Following  this  procedure  one  obtains  sufficient 
information  to  determine  if  the  observed  series  exhibits  stationary  or  non- 
stationary properties;  and  if  it  exhibits  non-stationary  properties  whether 
a first  or  second-order  difference  equation  will  remove  the  non-stationarit- 
ies. 

Once  we  have  obtained  a model  for  the  stationary  series,  a "backward 
filter"  is  applied  to  the  fitted  model  so  that  future  values  of  the  observed 
series  can  be  forecasted. 

2.2  Parametric  Time-Series  Models 

To  be  able  to  forecast  values  for  an  observed  series  we  fit  parametric 
time-series  models,  either  an  autoregressive,  a moving  average,  or  a combina- 
tion of  the  two.  These  stationary  stochastic  models  assume  the  process 
(series)  remains  in  equilibrium  about  a constant  mean  level.  The  general 
autoregressive  process  is  given  by: 

Xt  - U = ^(X^  - y)  + + c*,,, ( Xt_m  - p)  + Zt,  (2.7) 

where  p is  the  mean  of  X^,  is  a purely  random  process,  [14]  , and  m is  the 
order  of  the  process.  The  general  moving  average  process  is  given  by: 
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(2.8) 


Xt  - u = Zt  - BiZ^j 6qZt_q; 

u and  Zt  are  as  defined  above,  and  q is  the  order  of  the  process.  The 
general  mixed  autoregressive-moving  average  process  is  given  by: 
xt  - u = - y)  + •••  + o^n ( Xt_H,  - y)  + Zt 

" 6lZt-l 6qZt-q  * (2*9) 

where  q is  independent  of  m. 

We  shall  now  consider  the  criterion  for  selecting  the  process,  its  order 
(which  gives  the  best  fit  to  an  observed  series),  the  procedure  to  estimate 
its  parameters,  diagnostic  check  of  goodness-of-fit,  and  how  the  model  can  be 
employed  in  forecasting. 


To  estimate  the  parameters  for  the  autoregressive  process,  [l4j , we 
first  assume  that  the  Z*  process  is  normal  with  zero  mean  and  variance  a\. 
Then  the  log-likelihood  function  for  fixed  m,  conditional  on  the  values  xj, 
x2*  ***  * ^ can  expressed  as 


l(p,  aj,  **•  , | Xj,  •••  , Xjj,)  = - (n  - m)  (In  /(2n)  + In  a2) 

" 2^  t=L  f(Xt  ' U)  " ai(xt-l  " w) “m  (xt-m  - *• 


(2.10) 


For  estimating  the  parameters  y,  cij  , •••  , o^,  we  need  only  consider  the  sum 
of  squares  function 

S(y.  otj  , •••  , | Xj,  •••  , xm) 

= tJ+1  [(Xt  ' y)  " al(xt-l  - w)  - *•*  - " »*)]*.  (2*U) 

Now  assuming  that  y may  be  approximated  by  X and  that  the  sample  autocovari- 


ance function  (2.3) 


Cx*(j)  “ t=m+l  (Xt  ' n)(xt-J  " ^ j = *“  * m’ 

then  the  maximum  likelihood  equations  may  be  expressed  as 

cxx(j)  = SlCxx(j  - 1)  + S2cxx(j  - 2)  + ...  + - m). 


(2.12) 


(2.13) 


where  j = 1,  2,  .*•  , m.  Solving  the  m simultaneous  equations  one  obtains 
the  estimates  oq  , «••  , ctp,.  The  residual  sum  of  squares  may  be  expressed  as 
S(y,  alt  .*.  , am)  ~ (n  - m)[cxx(0)  - »1cxx(  1) Vxx(m)],  (2*14) 


and  the  residual  variance  by 


=1  ■ ST-k-n  sec.  s,  . - . %>. 


(2.15) 


To  estimate  the  parameters  for  the  moving  average  process  and  the  mixed 
autoregressive-moving  average  process,  we  use  a numerical  technique  to  build 
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up  the  log-likelihood  function  recursively,  [14],  By  varying  the  values  of 
the  parameters  (usually  between  -1  and  +1),  we  can  search  for  the  parameter 
estimates  which  minimize  the  sum  of  squares  function  for  each  process.  For 
example,  for  the  general  moving  average  process  the  sum  of  squares  function 
is  given  by: 

S(0,  Bx  . . 8q)  = z*  , (2.16) 

t-q 

where 

zt  = xt  - P + f^t-l  + •••  + Bqzt_q  , 
and  zt  = 0 for  t < q.  The  residual  sum  of  squares  is  given  by 

S|(q)  = S(y , Bj,  •••  , Bq)/(n  - q - 1).  (2.17) 

Similarly,  the  residual  sum  of  squares  for  the  mixed  autoregressive-moving 
average  process  can  be  expressed  as 

s|(m,  q)  = — - S(y,  cm  , , o^,  , •••  , fL).  (2.18) 

n-Zm-q-1  1 h 

2.2.3  Checking  the  Fit 

Once  we  have  fitted  a model  to  the  stationary  series,  we  must  determine 
the  adequacy  of  the  model.  If  it  was  found  necessary  to  filter  the  observed 
series,  the  first  step  is  to  apply  a "backward  filter"  of  the  same  form  so 
that  the  fitted  model  represents  the  observed  series.  Thus,  with  the  "back- 
ward filter"  inserted,  the  fitted  model  will  simulate  the  behavior  of  the 
observed  series.  Then  the  residuals— -the  observed  series  minus  the  modeled 
series— should  behave  approximately  like  a purely  random  process  (white 
noise),  that  is,  the  sample  autocorrelation  function  (2.3)  should  be  effec- 
tively zero  for  all  lags  except  the  zeroth.  To  determine  the  fit,  a test 
for  white  noise  is  applied,  [Mj. 


f 


T ■ 


I i 

' 
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2.2.4  Forecasting 

After  checking  the  fit,  we  can  use  the  resulting  equation  to  forecast 
a value  xt+^  , z >_  1,  when  we  are  currently  at  time  t.  This  forecast  is  said 
to  be  made  at  origin  t for  a lead  time  Z.  The  minimum  mean  square  error 
forecast,  [l3],  for  any  lead  time  Z is  given  by  the  conditional  expectation 
Et  xt+£  » xtU  at  origin  t,  given  knowledge  of  all  the  x's  up  to  time  t; 
that  is: 

xtU)  = Lxt+ll  * (2.19) 

The  required  conditional  expectation  occurring  in  the  forecasting  models  can 
be  found  using  Box  and  Jenkins  [l3j : 

Et[xt+jl  ” xt(j)»  ^t£zt+j]  = ® J = 1»  2,  (2.20) 

and 

Et[Xt-j]  = xt-j’  Et[zt-jl  = zt-j  j = °*  2*  '*•  (2-21) 

3.  IDENTIFICATION  OF  CLIMATOLOGICAL  DATA 

The  initial  step  in  the  analyses  of  the  two  observed  time  series  (daily 
rainfall  for  Greenwood  Lakes,  N.  J.,  and  Long  Branch,  N.  J.)  is  to  determine 
if  they  are  either  stationary  or  non-stationary.  Both  series  were  plotted 

in  an  attempt  to  graphically  detect  any  non-randomness  or  trend.  Figures  3.1 
and  3.2  display  the  365-day  rainfall  totals  for  both  stations,  and  certainly 
appear  to  exhibit  non-stationarities.  The  inference  of  the  graphic  displays 
were  statistically  tested  for  trend  using  Kendall's  Tau  test  as  detailed  in 
section  2.  For  further  evidence,  the  sample  autocorrelation  functions  were 
also  calculated  for  first,  second,  and  third  order  difference  filtered  data. 

The  critical  value  for  Kendall's  Tau  test,  flO]?  at  the  a = .05  level 
of  significance  is  + 1.645.  The  results  of  the  tests  are  given  in  table  3.1 
showing  a higher  order  filter  requirement  for  the  series  for  the  lesser 
rainfall  density  (Long  Branch,  N.  J.). 

* Chapter  5 
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FIGURE  3.2  MEASURED  DAILY  PRECIPITATION  FOR  LONG  BRANCH 

JANUARY  THROUGH  DECEMBER,  1974 


Table  3.1 

Kendall's  Tau 

Statistics  for  Trend,  Z 3 

0 

+1.645 

CALCULATED  STATISTICS 

Hq:  NO  TREND 

FILTER 

ORDER 

GREENWOOD 
LAKES  DATA 

LONG  BRANCH 
DATA 

GREENWOOD 

LAKES 

LONG 

BRANCH 

0 

9.158 

11.316 

Reject 

Reject 

1 

5.221 

7.321 

Reject 

Reject 

2 

1.532 

4.124 

Accept 

Reject 

3 

1.636 

Accept 

This  evidence  clearly  indicates  that  both  of  the  observed  series  exhibit 
non-stationary  properties  and  that  the  second  difference  data  for 
Greenwood  Lakes  showed  no  trend  at  the  a - .05  level.  Significiantly, 
however,  the  data  for  Long  Branch,  N.  J.,  required  a third  difference  filter 
to  remove  the  non-stationary  components. 

The  sample  autocorrelation  functions  of  both  the  observed  series  (see 
figures  3.3  and  3.4)  failed  to  dampen  out  rapidly.  This  indicator  further 
confirms  the  fact  that  the  data  is  non-stationary.  Figures  3.5  and  3.6  show 
the  sample  autocorrelation  of  the  first  difference  information  for  both 
series.  Here  again,  there  appears  to  be  no  rapid  dampening  to  indicate 
stationarity.  Graphic  displays  of  the  second  difference  data  for  both 
stations  are  shown  in  figures  3.7  and  3-8.  Here,  the  Greenwood  Lakes  rain- 
fall data  seems  to  exhibit  dampening,  but  admittedly,  only  slightly  more  than 
the  first  difference  data  of  figure  3.5.  The  Long  Branch  data,  on  the  other 
hand,  shows  little  tendency  toward  dampening.  Figure  3.9  shows  the  autocor- 
relation function  of  the  third  difference  information  for  Long  Branch,  N.  J. 
In  this  case,  dampening  becomes  strongly  evident  by  the  larger  peaks  in  the 
0-50  lag  range  and  smaller  peaks  in  the  270-365  lag  range.  These  particular 
characteristics  indicate  that  the  second  difference  series  for 
Greenwood  Lakes,  N.  J.,  and  the  third  difference  series  for  Long  Branch, 

N.  J.,  have  reached  statistical  equilibrium.  The  closeness  of  the  displays 
of  figures  3.5  and  3.7  for  Greenwood  Lakes  should  serve  as  an  indicator  that 
there  is  a strong  need  for  statistical  testing  (see  table  3.1)  along  with 
graphic  displays  of  the  data.  With  these  filtered  series,  one  can  now 
proceed  to  fit  forecasting  models  to  the  climatological  information. 
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FIGURE  3.3  SAMPLE  AUTOCORRELATION  OF  THE  DAILY  PRECIPITATION 
FOR  GREENWOOD  LAKES,  N.  J.,  1974 


FIGURE  3.5  SAMPLE  AUTOCORRELATION  OF  THE  FIRST  DIFFERENCE  DATA 

OF  THE  OAILY  PRECIPITATION  FOR  GREENWOOD  LAKES,  N.  J.,  1974 


FIGURE  3.7  SAMPLE  AUTOCORRELATION  OF  THE  SECOND  DIFFERENCE  DATA 
OF  THE  DAILY  PRECIPITATION  FOR  GREENWOOD  LAKES,  N.  J.,  1974 


E AUTOCORRELATION  OF  THE  THIRD  DIFFERENCE  DATA  OF  THE 
DAILY  PRECIPITATION  FOR  LONG  BRANCH,  N.  0.,  1974 


3.1  Fitting  the  Models 


To  fit  stationary  stochastic  models,  either  AR,  MA,  or  ARMA,  to  the 
filtered  information  as  outlined  in  section  2,  it  is  necessary  to  estimate 


the  parameters,  a-| , ag,  ...,  am>  6-j , 82*  . ...  8q,  for  each  process  and  for 


each  order  of  the  process  considered.  Following  the  proposed  procedure,  the 
parameters  were  estimated  with  the  restriction  that  they  lie  between  -1  and 
+1  to  insure  the  stationarity  and/or  invertibil ity  of  the  filtered  stochastic 
processes.  Simultaneously,  the  residual  sums  of  squares  were  computed  and 
divided  by  the  appropriate  degrees  of  freedom  to  obtain  the  residual  vari- 
ances. Figures  3-10  and  3.11  display  the  residual  variance  as  a function  of 
model  order  (m,q)  for  the  Greenwood  Lakes  data  and  the  Long  Branch  data, 
respectively.  The  interpretation  of  the  order  (m,q),  with  respect  to  the 
AR,  MA,  and  ARMA  processes,  is  discussed  in  section  2.2.  Clearly, 
in  both  cases  the  minimum  residual  variance  criterion  corresponds  to  mixed 
model  consisting  of  second  order  autoregressive  and  of  second  order  moving 
averages  components.  Thus,  the  order  (m,q)  for  both  the  Greenwood  Lakes  and 
the  Long  Branch  information  is  (2,2).  Using  the  corresponding  parameters, 
the  following  difference  equations  were  obtained  for  the  filtered  series: 

a.  for  the  Greenwood  Lakes  rainfall  data: 

(wt  - .0008)  = -0.582  (wt_1  - .0008)  - 0.227  (wfc_?  - .0008)  + 

+ 0. 762Zj._i  + 0.  172Zj._2»  aod,  (3.1) 

b.  *or  the  Long  Branch  rainfall  data: 

(ut  + .0022)  = -0.609  (ut  + .0022)  - 0.172  (ut  + .0022)  + Zt 

+ 0.769Zt_1  + 0.082Zt_2  (3.2) 

Since  a second  order  difference  equation  was  used  to  filter  the 
Greenwood  Lakes  data,  and  a third  order  difference  equation  was  used  to 
filter  the  Long  Branch  data,  it  is  now  necessary  to  make  use  of  the 
backwards  filters , as  outlined  in  section  2. 

The  filters,  respectively,  are: 


w„ 


and 


xt  " 2xt-l  + xt-2 


xt  " 2xt-l  + 2xt-2  ' xt-3 
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Order  - m,q 

FIGURE  3.10  MODEL  ORDER  vs.  RESIDUAL  VARIANCE  FOR  THE  DAILY  PRECIPITATION  1974 

FOR  GREENWOOD  LAKES.  N.  J. 


FIGURE  3.11  MODEL  ORDER  vs.  RESIDUAL  VARIANCE  FOR  THE  DAILY  PRECIPITATION,  1974 

FOR  LONG  BRANCH,  N.  J. 


We  insert  the  filters  into  equation  (3.1)  and  (3.2)  to  obtain  the 
appropriate  forecasting  models  that  will  be  used  to  characterize  the 
climatological  information.  Thus,  equations  (3.1)  and  (3.2)  become, 
respectively: 

xt  * 1.418xt_i  - 0.063xt_2  - 0.128xt_3  - 0.227xt_4  + .0015  + Z 

+ 0.762Zt_1  + 0. 1 72Zt  2 for  Greenwood  Lakes,  (3.3) 

and 

;t  * 2.391xt_1  - 1.345xt_2  - 0.311xt_3  + 0.093xt_4  + 0.172xt_s 

- .0039  + Zt  + 0.769Zt_^  + 0.082Zt_2  Lon9  Branch.  (3.4) 

Setting  the  unknown  Zt's  equal  to  their  conditional  expectations  of  zero  and 
assuming  the  values  x^,  xt_2>  ....  have  been  realized,  one  can  use 

equations  (3.3)  and  (3.4)  to  simulate  the  observed  climatological  series. 
In  addition,  if  t is  replaced  by  t + 1 in  the  above  equations,  one  can  fore- 
cast % steps  ahead,  l - 1,  2,  ....  L,  for  both  series.  Figures  3-12  and  3.13 
show  the  simulated  information  for  both  stations  which  clearly  fits  the 
observed  series  very  well. 

Tables  3-2  and  3.3,  p.  28  show  the  l step  ahead  forecasts  (up  to  l = 11 
lead  times  in  advance)  at  origin  t = 142  and  t = 144,  respectively,  with  the 
associated  confidence  intervals,  and  updating  of  the  forecasted  values. 
Ordinarily,  as  l increases,  the  forecasts  become  less  accurate.  The  short- 
term accuracy,  however,  can  be  maintained  by  updating  the  forecasted  values 
of  the  series  as  additional  information  becomes  available.  For  example,  the 
t = 142  origin  forecast  of  x^44  may  be  updated  to  become  the  t = 143  origin 
forecast  of  x^44  by  adding  a constant  multiple  of  the  one-step  ahead  forecast 
error,  e^Z^  3 0-]  Z^,  to  the  t 3 142  origin  forecast  of  x-|44-  The  fore- 
cast error  for  this  case  is, 

Z1 43  * X1 43  ' x143  * 

•and  6Z  = 63  is  given  by  ® 4^-83.  This  is  done  when  xt+1  = x143  becomes 
available.  The  basis  for  updating  the  original  forecasted  values  for  l steps 
ahead  as  additional  observations  become  available  is: 

xt+l^)  = xt^  + ^+%^t+l  * (3.5) 
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SIMULATED  PRECIPITATION  DENSITY  SERIES  USING  THE  MIXED  (ARMA) 
'EL  vs.  THE  OBSERVED  SERIES  FOR  GREENWOOD  LAKES,  N.  J.,  1974 


00 

o 


ll^tuey  j.o  saqoui  }uaiBAinb3 


Table  3.2  Forecasted  Values  of  Rainfall  Density  for 
Greenwood  Lakes,  NO  (1974)  at  Origin  t=142  and  Updating 
Under  the  Assumption  That  x143  Becomes  Available 


DATE 

ACTUAL  VALUE 
( Inches) 

LEAD 

TIME 

FORECAST 
( Inches) 

95%  PROBABILITY 
LIMITS 

UPDATED  FORE- 
CAST (Inches) 

5/22/74 

0.02 

-- 

— 

5/23/74 

0.15 

1 

0.148 

£ .327 

— 

5/24/74 

0.35 

2 

0.344 

£ .569 

.347 

5/25/74 

0.00 

3 

0.000 

± .832 

.005 

5/26/74 

0.00 

4 

0.007 

£ 1.128 

.001 

5/27/74 

0.05 

5 

0.063 

£ 1.451 

.055 

5/28/74 

0.06 

6 

0.050 

£ 1.800 

.060 

5/29/74 

0.13 

7 

0.127 

£ 2.172 

.138 

5/30/74 

0.00 

8 

0.000 

£ 2.566 

.013 

5/31/74 

1.10 

9 

1.101 

£ 2.566 

1.115 

6/01/74 

0.05 

10 

0.020 

£ 2.566 

.036 

6/02/74 

0.00 

11 

0.000 

£ 2.566 

.018 

Table  3.3  Forecasted  Values  of  Rainfall  Density  for 
Long  Branch,  NJ  (1974)  at  Origin  t=144  and  Updating 
Under  the  Assumption  That  x^45  Becomes  Available 


DATE 

ACTUAL  VALUE 
( Inches) 

LEAD 

TIME 

FORECAST 
( Inches) 

95%  PROBABILITY 
LIMITS 

UPDATED  FORE 
CAST  (Inches 

5/24/74 

0.02 

- - 

5/25/74 

0.00 

1 

0.007 

1 .920 

5/26/74 

0.02 

2 

0.018 

i 1.974 

.030 

5/27/74 

0.11 

3 

0.112 

± 3.516 

.156 

5/28/74 

0.01 

4 

0.011 

* 5.602 

.009 

5/29/74 

0.00 

5 

0.010 

* 5.602 

.000 

5/30/74 

0.00 

6 

0.002 

* 5.602 

.001 

5/31/74 

0.25 

7 

0.251 

* 5.602 

.253 

6/01/74 

0.58 

8 

0.575 

* 5.602 

.595 

6/02/74 

0.00 

9 

0.000 

± 5.602 

.000 

6/03/74 

0.00 

10 

0.037 

* 5.602 

.012 

6/04/74 

0.00 

11 

0.000 

* 5.602 

.000 

28 
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To  check  how 

the  residuals  were  calculated  using  z^  55  xt.-  xt  . The  simulated  (one-step 
ahead  forecast)  xt  is  subtracted  from  the  original  series,  xt.  Next,  the 
sample  autocorrelation  function  of  the  resiudals,  for  lags  of  0-364,  was 
computed  according  to  equations  (2.3)  and  (2.4)  . For  the  models  to  fit 
the  observed  information  well,  the  sample  autocorrelation  function  should  be 
effectively  zero  for  all  but  the  zero1"*1  lag.  When  the  observed  series  is 
sufficiently  large  (n  > 50),  the  sample  autocorrelation  of  the  residuals 
rzz(k)  ^ N(0,  Vn),  [11].  Thus,  for  both  stations,  the  standard  deviation 
of  the  sample  autocorrelation  is: 

— - = 0.0523 
/ 365 

and  the  95%  confidence  limits  are: 

+1 .96  (0.0523)  = +0.1025  . 


One  would  expect  that,  at  the  5%  level  of  significance,  365(.05)  or  19  of  the 
sample  autocorrelations  will  lie  outside  of  the  above  limits.  Hence,  from 
the  results  of  the  sample  autocorrelation  of  the  residuals,  one  can  conclude 
that  the  models  fitted  to  the  Greenwood  Lakes  and  Long  Branch  climatological 
data,  equations  (3.3)  and  (3.4)  , give  a good  representation  of  the  values 
realized  for  1974.  Tables  3.4  and  3.5  verify  these  conclusions. 


3.3  Proposed  Approach  vs.  Aka ike's  Approach 


Another  approach  to  the  classification  of  time-series  models  is  intro- 
duced by  Akaike,  [12],  which  uses  the  criterion  of  final  prediction  error 
(FPE).  This  method  was  used,  [8],  [6],  in  the  determination  of  "best" 
models  for  climatological  information.  Akaike  assumed  that  an  AR  model 
(equation  2.7)  has  a zero  mean  and  observations  acquired  from  points 
equally  spaced  in  time.  Using  the  Yule-Walker  equations  with  the  estimate 


t > m the  mean  square  of  residuals 
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TABLE  3.4  Sample  Autocorrelation  of  the  Residuals,  rzz^,  for  the 
Simulated  Precipitation  Density,  Greenwood  Lakes,  N.  J.  (1974) 
Confidence  Interval  = +0.1025 


Lag  K 


r — nrr 

Sample  Autocorrelation,  zzv  ' 


1-10 

1 .000 

r 1 49 

r408 

.016 

r002 

.039 

r064 

.108 

.037 

r061 

11-20 

r 009 

.009 

.052 

rllO 

r048 

.179 

r054 

r046 

.009 

.029 

21-30 

.063 

t053 

rOll 

r026 

r050 

.056 

.008 

.076 

r003 

T 1 00 

31-40 

.073 

rOl  4 

.003 

r080 

r030 

.182 

r036 

r045 

r032 

.054 

41-50 

r009 

t1  31 

.139 

.039 

r073 

.008 

r026 

.081 

rOl  3 

r031 

51-60 

.045 

r095 

r061 

.116 

.095 

r061 

r087 

.043 

.017 

r031 

61-70 

.019 

t025 

.016 

.026 

.025 

r004 

r061 

.075 

r045 

r 1 24 

71-80 

.097 

.086 

.009 

r 1 06 

.050 

.070 

r 1 30 

r055 

.074 

.127 

81-90 

r028 

t!02 

.054 

rOOl 

r072 

.031 

.083 

r029 

T 1 01 

.056 

91-100 

.042 

r027 

.006 

.032 

r 1 05 

.076 

.028 

r039 

.059 

r066 

101-110 

t024 

.057 

.013 

r048 

r023 

r036 

.097 

.039 

r053 

rOl  1 

111-120 

.011 

r044 

rOl  8 

.070 

.002 

r035 

.019 

.035 

.047 

r089 

121-130 

.032 

.147 

.012 

r092 

r052 

.020 

.026 

.009 

.020 

r039 

131-140 

.014 

r 004 

.006 

.011 

r033 

r022 

.004 

.079 

.009 

.083 

141-150 

r 005 

.041 

.007 

r028 

.002 

.027 

.008 

r057 

r024 

.098 

151-160 

rOll 

t045 

.013 

r002 

r024 

rOOl 

.057 

.001 

r028 

r025 

161-170 

.007 

.036 

.026 

r024 

.028 

.045 

r063 

r025 

.066 

.002 

171-180 

rOll 

tOI  5 

r006 

.024 

.007 

.063 

.029 

r 041 

.027 

.009 

181-190 

t003 

r01  4 

r044 

.017 

.074 

.013 

r040 

.011 

.014 

r058 

191-200 

t003 

.096 

r024 

r035 

.007 

.004 

.001 

r032 

.004 

.046 

201-210 

t004 

t021 

.001 

.005 

.007 

.011 

.009 

.008 

r002 

r003 

211-220 

t004 

.004 

.002 

r003 

r002 

rOOl 

.001 

rOOl 

.001 

.001 

221-230 

rOOl 

rOOl 

.001 

rOOl 

.0003 

.0004 

.0004 

.0003 

.0003 

.0002 

231-240 

.0001 

.0009 

r0003 

r0003 

r0003 

rOOOl 

r0002 

r0003 

r0003 

r0003 

241-250 

.0002 

rOOOl 

rOOOl 

r0002 

.0001 

.0002 

rOOOl 

.0009 

r0007 

r0009 

251-260 

.0001 

r0008 

r0004 

r0003 

r0003 

r0004 

.0002 

r0002 

r0004 

r0003 

261-270 

r0002 

t0003 

r0003 

r0002 

rOOOl 

r0002 

r0004 

r0003 

r0003 

r0003 

271-280 

r0002 

.0002 

.0002 

r0005 

r0003 

r0002 

r0003 

r0002 

r0002 

r0002 

281-290 

r0002 

r0002 

r0002 

r0002 

r0002 

r0002 

r0002 

r0002 

.0001 

.0003 

291-300 

tOOOI 

r0004 

r0002 

r0002 

r0002 

r0002 

r0002 

r0002 

rOOOl 

rOOOl 

301-310 

t0002 

rOOOl 

rOOOl 

rOOOl 

rOOOl 

r0002 

rOOOl 

rOOOl 

rOOOl 

rOOOl 

311-320 

.0001 

.0002 

.0001 

.0001 

.0001 

.0001 

.0002 

.0001 

.0004 

.0001 

321 -330 

t0002 

rOOOl 

.0001 

.0001 

.0000 

.0000 

rOOOl 

rOOOl 

rOOOl 

rOOOl 

331-340 

.0000 

rOOOl 

rOOOl 

.0001 

.0001 

.0003 

.0004 

.0002 

r0002 

rOOOl 

341-350 

tOOOI 

.0000 

.0005 

.0002 

r0005 

r0002 

rOOOl 

.0001 

.0001 

0001 
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TABLE  3.5  Sample  Autocorrelation  of  the  Residuals,  rzz^,  for  the 
Simulated  Precipitation  Density,  Long  Branch,  N.  J.  (1974) 
Confidence  Interval  = +0.1025 

Lag  K Sample  Autocorrelation,  r2Z(k) 


1-10 

1.000 

r249 

,246 

.156 

.036 

.055 

.115 

,046 

.090 

.074 

n-20 

.034 

.039 

.086 

r097 

.125 

.149 

r028 

.015 

.039 

.069 

21-30 

.061 

rOOl 

.077 

.019 

.049 

.031 

.061 

.048 

,007 

.074 

31-40 

.029 

.058 

.037 

,01  9 

.087 

.046 

.015 

.053 

,002 

.070 

41-50 

.039 

.003 

.054 

.005 

.099 

.005 

r019 

.076 

.050 

.013 

51-60 

.027 

.034 

.011 

.045 

.058 

tOI  7 

.068 

.001 

.040 

.073 

61-70 

t145 

.186 

.049 

,040 

.034 

.017 

.021 

.057 

.014 

,008 

71-80 

.039 

.065 

,044 

.067 

.019 

.009 

.030 

,029 

.067 

.050 

81-90 

,01  6 

.004 

.039 

.041 

,063 

.060 

.086 

,053 

.015 

.047 

91-100 

,032 

.067 

rOl  1 

.004 

.027 

.052 

,064 

.070 

.038 

,114 

101-110 

.088 

.087 

r065 

.027 

,073 

.108 

.041 

,029 

,106 

.126 

111-120 

.049 

t 071 

.025 

.031 

r 1 1 6 

.142 

,003 

,001 

,033 

,024 

121-130 

.080 

r023 

rOl  2 

.024 

,01  7 

.0003 

.034 

,049 

.041 

.010 

131-140 

,076 

.098 

r020 

r 1 01 

.090 

.010 

,027 

,006 

,013 

.002 

141-150 

.002 

tOII 

.015 

r019 

r039 

.043 

.004 

,026 

,015 

.004 

151-160 

.005 

rOll 

t008 

,005 

.003 

rOl  5 

.0008 

.007 

,022 

,005 

161-170 

.016 

r029 

.006 

,012 

.006 

.001 

,008 

,037 

.037 

,017 

171-180 

,032 

.039 

,01  6 

r034 

.031 

r054 

.043 

,0005 

,023 

,008 

181-190 

.006 

r017 

r025 

.036 

r005 

r023 

,008 

.002 

,007 

.001 

191-200 

,007 

,010 

,002 

,006 

.004 

r002 

,043 

.016 

.019 

,013 

201-210 

tOI  7 

rOll 

.004 

,0005 

r003 

tOI  2 

,011 

.006 

,002 

,011 

211-220 

,009 

,006 

,005 

,005 

r007 

t004 

,005 

,007 

,004 

,005 

221-230 

,006 

r005 

,005 

,006 

r006 

t006 

,005 

,005 

,006 

,005 

231-240 

,005 

,005 

,005 

r007 

r003 

r 004 

,005 

,004 

,005 

,006 

241-250 

,004 

,004 

,005 

,007 

r003 

r005 

,004 

,004 

,006 

,003 

251-260 

,004 

t005 

,004 

t005 

r004 

r005 

,004 

,004 

,004 

,004 

261-270 

t004 

,004 

r004 

,005 

r003 

r003 

,004 

,004 

,004 

,004 

271-280 

,005 

r003 

r003 

t004 

r004 

t004 

,004 

,004 

,004 

,003 

281-290 

t004 

r003 

,003 

t003 

-004 

t003 

,003 

,006 

,001 

,002 

291-300 

t004 

r003 

r003 

t003 

r003 

t003 

,003 

,003 

,003 

,003 

301-310 

r003 

r003 

t003 

,003 

r003 

t003 

,002 

,002 

,003 

,002 

311-320 

r002 

r002 

t002 

r002 

r003 

t002 

,002 

,002 

,002 

,002 

321-330 

r002 

,002 

t002 

r002 

r002 

r002 

,002 

,002 

,002 

,002 

331-340 

t002 

,002 

tO  02 

r002 

t004 

.0002 

.0006 

.002 

.001 

.002 

341-350 

r002 

,0005 

,0009 

r002 

rOOl 

rOOl 

,001 

,001 

,009 

,0004 

(3.6) 


1 N P 

Ro  = N Z (xt  " 1 ounxt-m) 

p N t»l  r m=l  t m 

are  to  be  minimized. 

The  FPE  criterion  depends  upon  the  use  of  the  mean  square  error: 

Sp  = N^P  Rp  (3.7) 

P+1 

and  (FPE)p  = 1 + — jq — , Sp  -v  Xp  • The  smallest  FPE  will  give  the  best  estima- 
tion of  the  parameters.  Table  3.6  below  shows  the  best  parameter  order  for 
both  approaches. 


Table  3.6  Best  Parameter  Order  for  the  Minimum 
Residual  Variance  and  Akaike's  Criteria 


MODEL 

ORDER  (m,q) 

RESIDUAL  VARIANCE 

Greenwood 

Proposed 

(2,2) 

0.163 

Lakes,  NJ 

Aka  ike (FPE) 

(8,0) 

0.158 

Long  Branch, 

Proposed 

(2,2) 

0.717 

NJ 

Akaike(FPE) 

(5,0) 

0.679 

The  difference  between  methods  is  obvious  from  the  above  table.  Based  on  the 
principle  of  parsimony,  [13],  one  should  reject  the  Akaike  classifications 
as  having  a substantially  larger  number  of  parameters  with  which  to  reckon. 

It  is  noteworthy  to  mention  that  in  Akaike's  method,  the  unbiased  estimate 
of  the  mean  square  error  will  usually  produce  good  predictions,  but  the 
variance  will  become  quite  large  in  the  analysis  of  the  spectrum.  Further, 
this  same  method  is  confined  only  to  the  autoregressive  models  and  does  not 
address  the  utility  of  the  moving  averages  and  mixed  models  in  analyzing 
climatological  information.  By  virtue  of  the  larger  number  of  parameters  of 
the  Akaike  classification,  the  number  of  previous  observations  required  to 
begin  accurate  forecasting  is  substantially  less  with  the  minimum  residual 
variance  classification.  For  instance,  in  the  case  of  the  Greenwood  Lakes 
information,  Akaike's  method  would  require  ten  previous  observations  for 
forecasting,  while  the  minimum  residual  variance  classification  would  require 
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only  four  previous  values.  The  fact  that  fewer  previous  observations  are 
required,  implies  that  a more  practical  near-real -time  forecasting  scheme  is 
possible  from  the  computational  point-of-view. 


| 

j 


3,4  Summary  and  Conclusions 

In  this  section,  the  procedural  approach  developed  in  section  2 and 
exercised  in  section  3 was  used  to  characterize  1974  climatological  informa- 
tion for  New  Jersey.  Specifically,  time-dependent  rainfall  data  for 
Greenwood  Lakes,  N.J.,  and  Long  Branch,  N.J.,  acquired  from  the  NOAA  was 
modeled  and  analyzed.  The  information  consisted  of  daily  rainfall  accumula- 
tion for  the  two  sites  taken  from  recording  rain  gauges  during  the  calendar 
year  1974.  The  two  sites  chosen  were  representative  of  northern  and  central 
New  Jersey  climate.  In  general,  the  climatological  characteristics  were 
fairly  uniform  within  a 15  Km  radius  of  each  station.  Therefore,  the  proce- 
dural approach  in  analyzing  this  type  of  data  would  be  relevant  for  use  by 
the  Army  intelligence  community  in  a tactical  situation. 

The  climatological  data  were  shown  to  be  non-stationary  realizations, 
and  following  the  procedural  approach  recommended  in  section  2,  the  following 
stochastic  processes  were  formulated  as  the  most  appropriate  characteriza- 
tions: 

i.  for  Greenwood  Lakes,  N.J.: 

xt  = 1.418xt_1  - 0.063xt_2  - 0.128xt_3  - 0.227xt_4  + 0.0015  + 1 
+ 0.762Zt_1  + 0.172Zt_2 
and 

ii.  for  Long  Branch,  N.J.: 

A 

xt  = 2.391xt_i  - 1.345x^_2  - 0.311xt_3  + 0.093xt_4  + 0.172x^  ^ 

- 0.0039  + Zt  + 0.769Zt_1  + 0.082Zt_2  . 

These  models  were  selected  on  the  basis  of  the  critierion  of  minimum  residual 
variance.  The  results  of  the  diagnostic  check,  through  simulating  the 
observed  series,  show  this  to  be  appropriate  with  respect  to  identifying  the 
actual  difference  equations  which  characterize  the  climatological  data.  As 
expected,  the  analysis  yielded  similar  difference  equations  for  both  sites 
(see  equations  3.1  and  3.2),  Furthermore,  we  have  structured  tables  that 
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show  short  and  long-term  forecasts  for  both  Greenwood  Lakes  and  Long  Branch, 
N.J.  In  these  cases,  ARMA  (2,2)  models  adequately  characterize  the  under- 
lying process  as  illustrated  by  figures  3-12  and  3-13,  and  in  tables  3.2  and 
3.3. 

The  information  gained  from  the  modeling  and  analysis  will  provide 
communicators  and  communications  planners  with  a mechanism  to  determine 

1 

short-term  future  comnuni cations  outages  in  the  GHz  range.  This,  in  turn, 
will  enable  planners  to  incorporate  suitable  alternatives  in  a tactical 
situation.  Systems  designers,  on  the  other  hand,  can  use  the  proposed 
modeling  technique  to  simulate  "worst  case"  outages  due  to  heavy  precipita- 
tion and  thereby  determine  adequate  equipment  operational  margins. 
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